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ABSTRACT 


NASA/Langley Research Center has recently embarked on a comprehensive 
program of evaluating higher harmonic bl ide feathering for helicopter vibra- 
tion reduction. Recent wind tunnel tests at the Langley Transonic Dynamics 
Wind Tunnel (TDT) have confirmed the effectiveness of higher harmonic con- 
trol in reducing articulated rotor vibratory hub loads. Several predictive 
analyses developed in support of the NASA program have been shown to be 
capable of calculating single harmonic control nputs required to minimize 
a single 4P hub response. In addition, a multiple-input, multiple -output 
harmonic control predictive analysis has been developed in support of an 
upcoming wind tunnel investigation of this mode of higher harmonic < mtrol. 
.Ml techniques developed thus far obtain a solution by extracting empirical 
transfer functions from sampled data. Algorithm data sampling and pro- 
cessing requirements are minimal to encourage adaptive control system 
application of such techniques in a flight environment. This paper is a 
status report on contract work to date. 
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Introduction 


The reduction of rotor-induced helicopter fuselage aeroelastic response 
carries with it the obvious benefits of enhanced operating life for critical 
aircraft components and reduced vibration of the controls and attendant pilot 
fatigue. Close examination of the rotor aerodynamic operating environment 
reveals the following: 1) rotor lift and drag and, hence, wake vorticity are 

a function of airstrcam velocity; Z) each rotor blade continually passes over 
its own wake and the wakes of other blades: 3) azimuthally varying vortex 
strengths and local airflows give rise to shear forces at the hub which are 
transmitted through the rotor shaft into the fuselage. Thus, blade root 
shears at harmonic multiples of the number of blades are the combined 
result of blade airloads, hub motion and various types of blade motion rela- 
tive to the hub. As these shears transmit from the rotating to the non- 
rotating system of a 4-bladcd helicopter, they give rise to the following: 

1) 3P and 5P flapwise root shears result in 4P hub pitching and 
rolling moments. 

2) 4P flapwise root shears feed into 4P shaft axial forces. 

3) 3P and 5P chordwise root shears produce 4P hub fore-aft and 
lateral forces. 

4) 4P chordwise root shears give rise to 4P hub yawing moments. 

It can be shown for an n-bladed rotor system that, assuming all blades see 
the same loading at the same azimuth position, only nW frequency harmonics 
will be seen in the fixed system. 

Thus, by utilizing phase relationships between various sources of rotor 
loads and introducing third, fourth and fifth harmonic pitch variations, the 
corresponding airload harmonics can be influenced. These airload harmon- 
ics, in turn, would influence blade motion and the modified motions and 
airloads together would modulate blade root shears. This is the premise 
upon which higher harmonic control is based, as illustrated in Figure 1. 

For a four-bladed helicopter, higher harmonic pitch control as a means of 
vibration reduction requires superimposing 4P swashplate motion upon col- 
lective and cyclic control inputs. When external 4P control of the swashplate 
is present, third, fourth and fifth harmonic pitch variations may be obtained 
independently of the pilot's use of collective and first harmonic pitch. Fourth 
harmonic pitch can be achieved by oscillating the swashplate verucally as a 
fourth harmonic about the collective position set by the pilot. The technique 
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Figure 1. Helicopter Vibratory Loads Reduction 
by Higher Harmonic Control 


for imposing third and fifth harmonic blade feathering involves fourth 
harmonic tilting of the swashplatc. Consider the following blade pitch 
schedule 6 (4>) as a function of lateral swashplate tilt 4> 


e (4<) = ^ (4i) sin 


( 1 ) 


If swashplatc tilt is allowed to vary as a fourth harmonic with amplitude A, 
then 


0 (4;) = A cos sin 4^ 

= — (sin 5 4^ - sin 34i ) . 

If swashplatc tilt is taken about a longitudinal axis, 
then, 


( 2 ) 


6 (4«) > A cos 4«p COB 4< 


3 


J (cos 5*|i 


cot 3*i<) . 


(3) 


Similar results would be obtained if swashplate tilt was a Hint, phased fourth 
harmonic. It is clear, therefore, that any phase and magnitude combination 
of third and fifth harmonic pitch is possible by proper choice of amplitude 
and phase of fourth harmonic swashplate tilting about longitudinal and lateral 
axes . 

The Piaterial that follows describes an investigation of several techniques 
developed under NASA contract NAS 1-14552 to solve for required harmonic 
pitch amplitudes and phases. All algorithms developed thus far compute 
input signals required to null shaft oscillatory forces by extracting empirical 
transfer functions from sampled data. The various techniques are contrasted 
in terms of accuracy and their potential for application in a closed loop 
active control network. 


Review of Previous Research 


Emerging as the only flight test investigation of higher harmonic control for 
vibration reduction conducted to date, a 1962 test program described in 
Reference 1 employed an open-loop second harmonic blade feathering control 
system installed on a modified UH-IA. Provisions were made for in-flight 
adjustments of both amplitude and phasing of 2/ rev blade feathering. 

Program objectives were threefold: 1) reduce vertical fuselage vibrations: 

2) reduce oscillatory blade loads; and 3) delay the onset of retreating blade 
stall. Although the investigation of second harmonic control as a means of 
stall prevention was inconclusive, the program did demonstrate that a 
reduction in vertical vibration at the aircraft eg and a reduction in oscilla- 
tory blade bending moments can be achieved simultaneously through proper 
application of harmonic fe^ljhering. 

1967 analytical investigation of higher harmonic pitch control to eliminate 
the transmission of oscillatory vertical forces from the rotor is detailed in 
Reference 2. 
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First the required pitch angles were found which eliminated oscillator/ lift 
loadings at all radial stations by compensating for the azimuthal variation in 
blade tangential velocity and downwash noimniformities. The assumed con- 
trol system was such that radial and azimuthal variations in pitch angle 
inputs were possible. 

Subsequent portions of the study were directed at computing pitch angle 
inputs at the blade root necessary to null oscillatory root shears. For this 
case, the constraint of eliminating oscilh.tory lift at all radial stations was 
relaxed. 

Conclusions reached on the basis of analysis of a UH- 1 A-configured rotor 
were that radial schedules of higher harmonic pitch angles required to 
eliminate all oscillatory airloads were a function of flight condition, and 
would be difficult to accommodate such radial pitq^ variations in a flight 
application. The study also concluded that a majority of transmitted oscilla- 
tory vertical root shears can be eliminated by using only pitch angle inputs 
at second and fourth harmonics on a UH-IA. 

An extensio*' of the work performed in Reference Z is described in 
Reference 3 in which oscillatory drag or inplane forces are addressed in 
addition to vertical vibratory' loads. Elimination of oscillatory inplane loads 
was analytically achieved by superposing third and seventh blade feathering 
(second and sixth harmonic swashplatc lilting motion) upon normal feathering 
requirements for trim. Reference 3 chose also to estimate the effect of 
higher harmonic pitch control inputs on performance. In concluding that the 
majority of vertical and inplane root shears can be eliminated simultaneously 
with second, third, fourth and seventh harmonic root pitch control, Refer- 
ence 3 first demonstrated the feasibility of multiple-input, multiple -output 
higher harmonic control applications. In addition, results of Reference 3 
support the contention that higher harmonic pitch control has little effect on 
rotor performance, a fact later to be borne out in wind tunnel investigations. 

Although analytical treatments of higher harmonic control have established 
its feasibility, much of what is presently known about the technique has 
evolved through wind tunnel test programs. References 4, 5, 6, and 7 sum- 
marize wind tunnel investigations of harmonic control as well as the results 
generated from such programs. 

References 4 and 5 describe a 1973 investigation, based on experimental fre- 
quency response data, of vibration compensation by periodic variations of the 
primary controls. Vibratory forces and moments to be compensated on the 
2. 29-m (7. 5 ft) hingeless model were vertical force and pitching and rolling 
moments. Using transfer function methods in the frequency domain to cal- 
culate required control inputs as well as the resultant blade response, it was 
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found bladr pitch variations required for vibration alleviation vary as a 
function of advance ratio; less than I de^^rce for 0, 2 < p and 

3 degrees for p = 0. 85. The -IP compensation control inputs were effective 
in attenuating 3P and 5P flapw'ise bending moments while 2P flap moments 
were largely unaffected by 4P blade feathering. It is noted that References I 
and 5 report there was no instrumentation to measure vibratory hub pitching 
and rolling moments, and that these quantities were obtained from flapwise 
bending moments. Thus, the effect of inplane forces, vertical shear forces 
and blade torsion have been ignored. In addition, the effect of hlade feather- 
ing on rotor performance was not addressed. 

Reference 6 and its extension Reference 7 describes what has proven to be 
the most comprehensive and encouraging wind tunnel investigation of har- 
monic control to date. By imposing 4P swashplate vertical motion as well 
as pitching and rolling motion on a 3, 05-m (10. 0 ft) Mach-scaled hingeless 
rotor, it was observed that vibratory hub load reduction carries with it no 
penalty in terms of blade flap bending or rotor perlormance. While Refer- 
ence 6 describes performance benefits of 20 percent with 2P feathering. 
Reference 7 reported no such improvement with 4P feathering. 

Paramount in the conclusions of Reference 7 is the fact that simultaneous 
reduction of several vibratory hub load components is possible due to 
superposition of harmonic pitch effects, thus underscoring conclusions ^ 
Refe-ence 3. ' * ^ 

Although not discussed in previous work, high frequency blade pitching 
carries with it increased control loads by virtue of the tennis racket effect. 

In contrast, Multicyclic Flap Control, Reference 8, wherein high frequency 
elastic torsional deformations generated by a mid-span ser\n flap are 
employed to null hub vibratory loads, has no such penally. Inertia-induced 
control loads, in fact, are greatly reduced since only a flap is being driven, 
not the blade. In addition. Reference 9 reports mid-span pitching is far 
more efficient than full span pitch control in reducing first mode flap bending 
excitation. Inherent design drawbacks of servo-flap contr 'l are g"eater 
mechanical comple.xity in terms of rotor blade/ cont rol sysiem design 
integ ration. 

Although the above studies as well as others not mentioned here, were 
successful in demonstrating the potential of harmonic control to e.xpanct the 
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flifiht envelope of the hingeless rotor there exists b need for further study in 
several areas. Specifically, this paper intends to address the following: 

• The effect of -tF bwashplate vertical, pitch and roll motion on 
hub forces and moments of an articulated rotor model in the 
wind tunnel. 

• General characteristics of several algorithms developed to 
predict required swashplate inputs. 

• Considerations for applying harmonic control in an active feedback 
network. 

Wind Tunnel Testing 


The solution techniques explored thus far have been developed in support of 
Lang’cy Research Center's comprehensive program of evaluating active 
higher harmonic control for reduction of helicopter aeroclastic response. 

The program embodies the phases shown in Figure I, The higher harmonic 
control concept was evaluated recently using a 2. 74-m (9. 0 ft) diameter aero- 
elastically scaled rotor model in the NASA/Langlcy Transonic Dynamics Wind 
Tunnel (TDT). Characteristics of the model rotor are given in Figure 3, 

Table I and Figure 4. The test was conducted such that manual phase and 
amplitude sweeps of 4F collective, longitudinal cyclic, and lateral 
cyclic feathering were made to explore the frequency response of the 


ACTIVt CONTROLS FOR REDUCTION OF 
HELICOPTER AEROELASTIC RESPONSE 


• ESTABLISH HIGHER HARMONIC CONTROL EFFECTIVENESS 

- EXPERIMENTAL PROGRAM 

- 4 BLAOEO ARTICULATED ROTOR 

- SWASHPLATE EXCITATION 

- OPEN LOOP 

- SINGLE INPUTS AT FIRST, THEN MULTIPLE INPUTS 

• DETERMINE CONTROL LAWS FOR ACTIVE SYSTEM 

- ANALYTICAL PROGRAM 

• ASS'SS EFFECTIVENESS OF ACTIVE SYSTEM USING 
WIND TUNNEL MODEL TESTS 

• VALIDATE CONCEPT THROUGH FLIGHT TESTS 


Figure 2. NLRC Active Higher Harmonic 
Control Program 
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Figure 3. Blade Rotating Natural Frequency 
Characteristics 

TABLE I. SUMMARY OF ROTOR SYSTEM PROPERTIES FOR THE 
2. 7-METEK (9 FOOT) DIAMETER RESEARCH MODEL ROTOR SYSTEM 

Rotor Type 

Number of Blades 

Rotor Diameter 

Blade Chord 

Spar Chord 

Solidity 

Airfoil Section 

Theoretical Blade Twist 

Hinj;e Offset, Coincident 

Root Cutout 

Blade Elastic Axis 

Blade Pitch Axis 

Center of Gravity 


Fully Articulated 
4 

279 -_m (110 in. ) 

10. 77 -cm (4. 24 in. ) 

4. 72 -cm (1. 86 in. ) 

0. 0982 
NACA 0012 
0 degrees 
7. 62-cm (3. 0 in. ) 

31. 75-cm (12. 5 in. ) 

25% Chord Nominal 
25% Chord 

Adjustable (Outer Blade 22% 
thru 37%) 
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Fijjure 4. Aeroelastic Experimental Rotor System 


six hub forces and moments. Combined rotor / fuselage responses were 
measured by a six-component balance mounted below the pylon. Hub 
vibratory loads as measured by the balance were corrected to co^ipensato 
for n_.T-scale swashplatc inertial forces. Blade and pitch link loads were 
also monitored. 

Figures 5 and 6 depict a representative ^'ariation of 4P hub^ormai response 
with 4P collective feathering inputs. Azimuthal and radial coordinates of the 
Figure 4 phasor plot represent 4P hub normal response pha#e und magnitude, 
respectively. .Sample 4P collective inputs were made at a constant ampli- 
tude of 1/2 degree, the data points being annotated with input phase. 

Baseline hub response represents the ambient hub response to be compen- 
sated prior to any feathering inputs. Figure 6 illustrates an optimization 
technique by which a local minimum can be reached by varying 4P input 
phase, and further refined by modulating 4P input amplitude. 

Preliminary wind tunnel results concerning the effect of 4P feathering 
inputs on 4»lade loads are presented in Figures 7 through 12, and discussed 
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in greater detail in Reference 10, It was found that although flapwise and 
edgewise bending moments were fairly insensitive to 4P collective inputs at 
the "optirrised" condition of Figure 6 (0. ii degree amplitude, iO degree 
phase), root torsional moments and, hence, pitch link loads were aggravated 
by such inputs, as illustrated in Figure 9. Specifically, the 4P content of 
blade torsional moment and pitch link load was highly amplified as a result 
of feathering induced inertial loads. Figures 10, II, and 1^. 

As noted in Figure 3, however, the first elastic torsion mode of the blade 
was above lOP at the design operating speed. Thus, a more torsionfilly 
compliant rotor may not have undergone such dynamic amplifications 
in pitch link loading. The substantial increase in SP content of torsiene.i 
moment and pitch link load, as depicted in Figures |0 and 11, can be attrib* 
uted to inadvertent mixing of 4P collective signals with swashplate pitching 
and rolling niotion. 

By way of conclusion, wind tunnel testing of a dynamically scaled model 
rotor has established that higher harmonic blade feathering can minimize 
rotor-transmitted hub vibratory forces and moments. Although such feather- 
ing inputs may result in higher blade and control system loads, judicious 
blade design may provide relief from any loading penalties. 

Predictive Analysis Overview 

The effectiveness of higher harmonic control having been established 
through wind tunnel testing described in the previous section, a deterniinis - 
tic approach to solving for required feathering inputs is developed next. The 
desire to implement higher harmonic control in an adaptive control system 
has created a set of unique requirements for solution techniques. Algorithms 
to be employed must be efficient so as to permit a high solution refresh rate. 
In addition, the algorithm must require a minimum of sampled data, both in 
quantity and type, so *3 .not to impose burdensome data acquisition require- 
ments. Working within these requirements, there e.xist several control 
philosophies that can be accommodated: 

• A single swashplate degree of freedom (d. o. f. ) is used to 
control a single hub response. 

• A single swashplate d.o, f. is used to control an aggregate of 
hub responses. 


• Multiple swashplatc d.o.f, 'a are uaed to control a single hub 
response. 

• Multiple swashplate d.o.f. 's are used to control multiple hub 
responses. 

In addition, any one of these techniques may be implemented using multiple 
frequencies (i. e. , 4F. 8P, and 1 blade pitching ), 

Several predictive techniques have been developed to treat the single input, 
single output higher harmonic control philosophy. These techniques were 
developed in support of the 1977 NLRC TDT teat program, and, conse- 
quently, results have been generated from available data. Techniques 
dealing with multiple-input, multiple-output harmonic control have similarly 
been developed for a planned March 1978 TDT test program. The various 
techniques will be developed and contrasted in the sections that follow. 

Single Input, Single Output Predictive Analyses 

Three-Point Technique 

In developing predictive algorithms to calculate optimal swashplate inputs, 
advantage was taken of the almost linear relationship between 4P feathering 
inputs and 4P hub oscillatory forces and moments. The first such approach 
developed requires a priori knowledge of baseline vibra.ion levels and two 
samples of oscillatory output in response to known 4P feathering inputs. By 
using 4P collective swashplatc inputs to minimize 4P hub normal forces, the 
procedure can be outlined as follows: 


1) Define feathering inputs and response quantities as phasors 
having magnitude and phase as: 


ZBL 


^Kzbl 


'^'ZBL^ 


Phasor representing 
4P component of 
baseline hub normal 
f^rce, with amplitude 
|F'ZBL| and phase 
relative to an index 

*ZBL 
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) 



Fhasor representing 
4P component of hub 
normal force in response 
to an arbitrary 4P collec- 
tive input 

Fhasor represe'iting 
4P component of hub 
normal force in response 
to an arbitrary 4P col- 
lective input different 
from above 


Oqi = Phaser representing 

the first arbitrary 
4P collective pertur- 
j^tion with magnitude 

®01 I 

" ^l®02 ’ ^***02^ Phaser representing 

a second arbitrary 4P 
collective perturbation 

2) Transform phase and amplitude to sine and cosine magnitudes; 


F = F cos <li 

ZCBL ZBL ^ZBL 


F ~ F sin 4> 

ZSBL ZBL ’^ZBL 


®oc, ' 


G - 6 sin 

”oSl 01 ^01 


etc. 


IS 



3 ) 


Assume there exists a plane, Fifiure 13, that describes the 
relationship between IP hub normal force cosine magnitude 
and 4P swashplate collective sine and cosine magnitudes. 

Assume a similar plane exists for the suie magnitude of 4P 
hub normal force. In addition, establish the following limitations: 

a) The F,-_ and F planes arc not parallel to each other 

b) Neither of the two planes arc parallel to the plane of zero 

response = F^^ = 0 ) 

c) The locus of points representing the intersection of the 
planes is not parallel to the plane of zero response, 

4) Write equations for the and F^s planes in terms of two 

arbitrary coefficients and baseline magnitudes. 


A 6 BO (F - F ) = 0 

OC OS ' ZC zcbl' 


DO ^ EO * F - F ) = 0 
OC OS ZS zsbl' 


(4) 

(5) 


By substituting two frequency response samples, the coefficients 
in Equation (4) and Equation (5) may be determined. 

5) Referring to Figure 13, F^q and F 23 in Equations (4) and (5), 
respectively, may now be set to zero. This yields equations 
for two lines in the zero-response plane, represented by line 
segments A^ and CD in Figure 13. Simultaneous solution of 
the two equations yields point P, whose coordinates are the 
sine and cosine niagnitudes of the 4P collective input needed to 
null the sine and cosine magnitudes of hub response. 

A schedule of 4P collective inputs and resulting 4P hub normal force 
responses for a particular wind tunnel trim condition is presented in 
Table II. Test conduct was such that 4P collective amplitude was held con- 
stant while phase was swept manually in near 45 degree increments. 

Table UI presents results using data from Table II with the 3-point tech- 
nique. Results were calculated based on combinatorial permutations of two 
of the eight available data points to check solution consistency. A single 
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Figure 13. Planar Relationships Between 4P Collective 
Input and 4P Shaft Axial Force Sin and Cos Components 


solution was then obtained by neglecting those solutions with phase greater 
than one standard deviation away from the mean phase. A new mean phase 
and amplitude was then calculated from the reduced set of solutions. 

Results generated from test data were fairly consistent with three exceptions 
noted. Close examination revealed that when 4P perturbations were made 
180 degrees apart, as were the 3 cases of Table III, the planes that resulted 
were vertical and coplanar, thereby defining an infinite set of non-unique 
solutions. Rather than impose a sampling scheme to preclude 180 degree- 
apart sampling, it was deemed more appropriate to develop a technique that 
has no such constrain^.? 
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table II. SCHEDUI.E O. 4^ COLLECTIVE INPUTS AND HUB RESPONSES 

H = 0. 30 C j /a = 0. 075 

m,Her BUSe EealHe.In, I„. HeUc„,ie. 


Point 


Response of Interest 


213 

Baseline/Normal Force 

1 

214 

Collective/Normal Force 

215 

Collective/Normal Force 

216 

Collective/Normal Force 

217 

Collective/Normal Force 

218 

Collective/Normal Force 

219 

Collective/Normal Force 

220 

Collectif^/Normal Force 

221 

Collective/Normal Force 


4p Input 
Amp 

(deg) 

4P Input 
Phase 

(deg) 

4 P Output 
Ainp 

N (lbs) 

4 P Output 
Phase 

(deg) 

0. 00 

0. 00 

(25. 80) 
114. 8 

44. 00 

0. 50 

27. 00 

(35. 40) 
157. 5 

-138. 00 

0. 50 

73. 00 

(46. 60) 
207. 3 

-59. 00 

0. 49 

114. 00 

(65. 00) 
288. 1 

-26. 00 

0. 49 

166. 00 

' (73. 90) 

328. 7 

13. 00 

0. 49 

-151. 00 

(76.40) 
339. 8 

41. 00 

0.48 

-104. 00 

(70. 90) 
315.4 

75. 00 

0. 49 

-50. 00 

(57. 10) 
254. 0 

121. 00 

0.49 

- 1 

-9. 00 

(45. 00) 
200. 2 

165. 00 



TABLK m. SOLUTION INPUTS BASED ON 3-POINT TECHNIQUE 

(DATA FROM TABLE II) 


Three-Point Higher Harmonic Solution 


Data from 
Cases 



4P Input Amplitude 

4P Input Phase 

0. 21 1 

28. 100 

0. 210 

28. 132 

0. 210 

28. 224 

0. 197 

26. 131 

0. 21 1 

28. 374 

0. 210 

28. 220 

A 0.210 

^ 0.251 

28. 084 
26. 286 

0. 251 

25. 268 

0. 251 

32. 355 

0. 361 

91.018 

0. 24 1 

23. 078 

0. 223 

28. 452 

0. 250 

24. 54 3 

0. 251 

32. 710 

0. 24 1 

40. 101 

0. 218 
0. 208 
0. 251 

13. 641 

29.913 A 

33.437 " 

0. 244 

31. 596 

0. 264 

28. 002 

0. 102 

87. 186 

0. 250 

33. 180 

0. 252 

33. 091 

0. 256 

32. 975 

0. 231 

34. 196 

0. 220 

34. 396 

0. 208 

34. 652 


Linear Higher Harmonic Solution 
-JP Input Amplitude IP Input Phase 


0. ZiZ 


30. I 10 
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Nonlinear (Six-Point) Technique 


The second predictive approach to be investigated replaces the planes of the 
forme' approach with second>order surfaces defined by equations 6 and 7: 


A0 ^ + D6 ^ Ce 0 ♦ D0 E0 

OC OS OC OS OC OS 


F0 G0 ^ + H0 0 10 *■ J0 

OC OS OC OS OC OS 


F - F 
ZC ZCBL 

F - F 
ZS ZSBL 


( 6 ) 

(7) 


The shape functions above are similar to cubic polynomials used in finite 
clement plate analysis: the lower-order terms are retained to improve the 
approximation while higher-order terms are eliminated to reduce the number 
of samples required to define coefficients A through J. 

In addition to a baseline condition, five san^ple swashplate inputs and result- 
ing hub responses are required to uniquely define the two shape functions. 
Once the 10 coefficients are defined, a Newton-Raphson iterative scheme is 
used to solve the set of nonlinear algebraic equations for the required 
swashplate inputs. The scheme takes the following recursive form: 


whe re 


tolit, Hell - i=l N 


( 8 ) 


{e|. = 


OC 

*os 





^«os * ° 


+ C0^^ 
OC 

+ E 


""os * ' 


H0 

OC 

J 


^0 





and N is the number of iterations. The analysis is initiated (i « o) by using 
solution inputs derived from the linear 3>point technique or a reasonable 
guess. A satisfactory tolerance on convergence can usually be mot in 
4*5 iterations. Table IV presents results generated from Table II data using 
the nonlinear approach. GxllM agreement is seen between the linear and 
nonlinear algorithms, indicating hub forces are more linear than nonlinear 
with control inputs of small amplitude. 

The need for five data points as well as a baseline does not preclude active 
control applications of a nonlinear approach. Reference I 1 details a 
minicomputer-directed helicopter active control network where up to 
24 samples of rotor state were required each revolution. In the event a 
particular control application was unable to sustain the rather low 
6 samples/rev rate, the process could be retarded to provide a new solu- 
tion every several revolutions. The lowered solution rate would virtually 
be undetectable in the cabin environment. 

The paramount drawback to this type of analysis lies in its arduous data 
processing requirements. The inversion of 5 by 5 matrices coupled with an 
iterative solution process could erode control loop response. Thus, a third 
predictive analysis was investigated which has no inherent sampling con- 
straints nor exhaustive data processing requirements. 

Two-Point Techtiiquc 

The third single-input, single output algorithm to be investigated requires a 
baseline and only one sample frequency response to provide a solution. The 
technique, developed by Langley, is based on the following tacit assumption: 
if higher harmonic partial response is defined as that portion of hub 


TABLE IV. SOLUTION INPUT BASED ON 6-POINT NONUNEAR 
TECHNIQUE (DATA FROM TABLE II) 

Nonlinear Higher Harmonic Solution 
Using Newton-Raphson Iteration 

Number of Iterations = 3 


Data from Cases 

4P Input Amplitude 

4P Input Phase 

214, 215, 216, 217, 218 

0. 223 

29. 514 
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o 


• scillatory force due eolvly to 4F feathering inputs (i. e.^ total response 
ninus baseline response), then 4H sample input phase leads harmonic 
artial respunre phase by a constant amount. The validity of this assump* 
tion uiing Table 11 data is established in Table V. Referring to Figure 14, 
the algorithm can be summarized as follows H) ^ 

1) Using phasor notation, vectorially subtract ^e baseline 

hub response of interest from the perturbation hub response. 


F 3 F - F 
ZHH Z1 ZBl. • 

• • 


(9) 


o 


TABLE V. ASSUMPTION: SAMPLE INPUT PHASE 

LEADS lUGHER HARMONIC CONTROL PARTIAL 
RESPONSE PHAS* BY CONSTANT AMOUNT 


,0 


• 

4P Collective Input Phase, dejj 

4P HHC Partial Response Phase 

A 

27 

222 

196 

73 

® 275 

202 


O 


1 14 

310 

197 

166 

358 

193 

O 



209 

• / 

190 


• 


256 

90 

194 

310 • 

147 

197 


e 


351 O 

185 

194 


0 


o 
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Figure 14. Jwo-Point Higher Harmonic 
Solution Technique 
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3) 


4 ) 

O 


Rotate the higher harmonic partial respc^se phasor ulfi*til it 
opposes_the baseline phasor. If has magnitude and 

phase ^ZHH> step requires a rotation of 

magnitude 6, 


ZHH 


"ZDL ■ 


<P 


By virtue of the assumption that the diff«Q*ence between 
input phase and harmonic control partial response phase is 
constant fur a given flight conditioi\) the required cunti^ol 

input phase may be written 

O 


rol 


o 




= A0 


od 


O o 

^°ZHH ® ^^ZDL 


ISO)) 

O 


0 


Allowing that harivonic control partial response magnitiidt' 
|FzhhJ ** lirnear with control input amplitude, the required 
control input amplitude for nulling baseline response* is 


o 

o 

O 


(10) 
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Z3 






o 




ZBL 


ZHHl 


( 11 ) 


Table VI presents results generated by the twu-point approach using wind 
tunnel data from Table II, The level of agreement between this approach and 
the previous techniques indicates it is the most likely candidate for control 
applications, given its generality and simplicity. 

Thus, it has been established that by virtue of the almost linear relationship 
between feathering inputs and hub oscillatory forces, uevei.'*! technirjues 


TABLE VI. SOLUTION INPUTS .:.\SED ON TWO-POINT TECHNIQUE 

(DATA FROM TABLE U) 


Two -Point Higher Harmonic Solution 


Data from Case 

4P Input Amplituda 

4P Input Phase 

214 

0. 210 

28. 156 

215 

0. 221 

21. 627 

216 

0. 206 

27. i43 




217 

0, 236 

31. 391 

218 

0. 249 

33. 527 

219 

0. 244 

29. 763 

220 

0. 221 

2o. 891 

221 

0. 202 

29. 222 


Linear Higher Harmonic Solution 


4P Input Amplitude 4P Input Phase 

0. 220 28. 79S 
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fxisl l«»r I mj{ -IP ronlrol iiipiil t«» niinimi/t* b.iHflitu- vibrations on ihv 

baKiH ot fiampi (i fiai.i. Ibc iiotcnt i<i I «tf liirsv tfchnicpit* s, as well as others 
to (>e applied in u multiple-input, inultipli' -output mode, will be discussed 
in th«- next f»»-ctn>n. 

Multiple Input, Multiple Output Predictive Analyses 

A primary consi<l»Tation in apprnat hin){ hi^jher harmonic control in a multiple- 
input multiple-output mode is that there exist only three independent swash- 
plate deijrees of freedom to minimize six hub vibratory responses. One 
approach is to minimize just the three hub responses that largely contribute 
to vertical fuselage response, namely, 4P verti*al forces, 4P fore-aft forces 
and 4P pitching moments at the hub. Alternatively, one could address the 
responses contributing to lateral fuselage response: IP side forces and 4P 

roliins and vawing moments. However many hub forces will be minimized, 
once the control objective has bi'en established the predi:tive algorithm must 
be developed. 

In developing a predictive technique for this mode of harmonic control, it is 
desirable to again ininit.iize sampling and data processing requirementf . 

Thus, to sli plv extend the 6-point nonlinear technique to 3-inputs and 
3-outputs, would require 15 input perturbations, and 48 FFT's or spectral 
analyses, thereby, proving too burdensome for adaptive control systems. 
Similarly, tb three-point technique would require 6 input perturbations and 
i\ FFT's in a 3-input, i-output mode. Even extending the two-point technique 
to such a mode would require 3 perturbations and 12 FFT's per solution 
ipdate. A linear regression analysis in the frequency domain, discussed 
in the next section, requires three perturbations and only six FFT's and as 
sue!'., hold much promise in multiple - input, multiple output applications. 

Linear Regression Predictive Analysis 


The analysis of a linear system with p inputs and a single output shall be 
considered first. By virtu,i of systcivi linearity, response in the time domain 
can be written as a sum of p individual component responses 


P 

Y(t) - ^ Y (t) 
k=l 


( 12 ) 
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For thr caar of a lin^’lc hul) output, say normal force, and input from 3 
swaahplate d.o.f^s Fq. 12 I'ecorncs 


F (t) = F (t) f F (t) > F (l) + F (t) (13) 

z z , z , z. z 

1 2 3 o 

where Fj.^ represents a baseline normal force response and F^j^(t) repre- 
sents normal force response to the input when other inputs are zero. 

The individual component responses, can he represented as the con- 

volution of the inputs X|^(t) with the impulse response functions h|^(t). Hence 
the component outputs arc >{iven as follows: 




/ 




r)dT 


(14) 


The Fourier transform of Eq. 14 yields 


F ^(f) = 
zk 


H. 


(0 


Fz 


X. (f) 
k 


(15) 


where F 2 j^(f) and Xj^(f) are the complex Fourio- transforms of Fjjj^(t) and x(t) 
respectively. Thus, the Fourier transfoim Fz (f) for the total output is, 

• X 

F (f) = H, (f)X ,(f) f H, (f) X (f) + H, (f) X Jf) f F (f) (16) 

z Ip 1 2p 2 3p 3 zo 

z z z 


or, 

P 

Hi. (OX.(f) F (f) (17) 

z k„ k zo 

k= i r 

z 



where is the complex Fourier transform of the baseline normal 

force response. 

Now if p independent observations of the inputs and outputs are available, 
then H|^p (f) k 3 1 ^ , p can be obtained by solving the linear simul- 

taneous equation 


F 

z 



■ 

k=l 




+ F (f) 
zo 


(18) 


(v = 1, 2, ... p) 


where v is the observation index. As outlined in Reference 13, for the 
case when only a single observation is available, as would be desirable 
for applications to higher harmonic feathering, we can select a set of p 
frequencies for which the input vectors. 


x,(f,) 


x,(f ) = 

1 V 


X,(f ) 
I P 


I 

) 


x,(f ) = 

2 V 


X2(f2> 


X,(f ) 

2 p 


. . . . X (f ) 
P V 


are linearly independent. We can then obtain an approximation 

H)(p. (f) by solving the simultaneous linear equations 
*’z 


flkp^(0 to 


P 

w = <‘'W '■’> 

k= 1 F 

z 


(v = 1,2 p) 
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Acknowk‘d(<in){ the existence. of additive noise n(t), the observed output F (t) 
is represented as 


F (t) = F (t) f n(t) 
z z 


UO) 


Taking Fourier transforms of equation iO, 


F (f) = 
z 


P 

'i H (0 X (0 ► n(0 
k=l "f 


( 21 ) 


and thus if we solve the following equation for the unknowns (f) (k = 

1, i, ... p) in place of the equation for H. (f) * 

•"f 

z 


(0 X. (f ) 

z V k,, k V 

k= 1 r 

z 


( 22 ) 


(v = 1 , 2, . . , p) 


we get an approximation to (f) contaminated with the noise n(f ). 

If we assume tnere exists a set of frequencies |fn{ which H|jp^(f£) = 
Hkp . k = 2, . . , p, f = 1, 2, . . , N, then 


= 1 H » n(f,) 

k= 1 F 


( 23 ) 


(2= 1, 2. ... N) 
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linplifii a problem of catimatinK multiple regression coefficients. 

Now consider the complex random noise variables n^fv = 1, L, , N), which 
are mutually independently distributed with first and second moments 


E(n^) = 0 
E(Ren = E(Im n 

V V 



(24) 


and the observations are made on variables , X , (v 

I, 2, .. , N) satisfying the relation v v Pv 


F (f ) 

z V 




H. 


k=l 


kv 


(25) 


(v = 1. 2. ... N) 


where the inputs Xi , X> .... X_ are considered to be observable non- 
. , ‘v Pv 

random variables. 

The probability density function L for the noise n^ is given as follows: 


-r\ I ' y 

- — 2 ~2 ^ 

\2tt(t / \ 2(t v=l 


(26) 


Substituting for n^, from equation 25, we obtain 




cxpl- 


2(r v^l 


P 

F (f ) - "S H. X (f ) 
z V ^ k k V 
k= 1 r 

z 


(27) 


29 


The maximum probability estimate, according to Reference 13, . 

flpj,- ) of (Hip . , Hp., ) is obtained by finding the 

values of H|^p which minimize t^c sum of Squares in the second paren- 
thesis of the right hand side of equation il. These values are given as 
follows: 


- - 


* 



(F . X,) 
z 1 


F 

z 


(F , XJ 
z 2 

F 

7 . 

• 

= d"‘ 


• 

H 


(F- . X ) 

Pf 


z p 

z 


_ — 


(2H) 


where 


(Xi, Xj) (X^, Xj) . 
(X,, X^MX^, x^) . 


(X . X 1 
p 1 

(X , X J 
p ^ 


(X,, X ) (X,, X ) . 

I p 2 p 


(X . X ) 
p p -• 


(29) 


and the inner products are defined by 


N 

(X., X.) = ^ X. X- (30) 

J ‘ IV gv 

v= 1 

(jl 2 ~ If 2, . . , p ) 

( )’!' denotes complex conjugate 
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and it is assumed that N S p and D ^ exists, and for out case, (p a 3), 
* FKT of swashplate collective inputs at frequency f^ 

3 FFT of swashplate lateral cyclic inputs at frequency f^ 

X^ * FFT of swashplate longitudinal cyclic inputs at frequency f^ 


If the above analysis is performed twice again, one each for hub pitching 
moment and axial force, we obtain 



(M , X,) 

_y ‘ 

(M , X,) 

y 2 


. (M , X ) , 

L y p j 


H. 


-1 


(F . X,) 
X 1 

(F . XJ 

X I 


(F , X ) 
. X P . 


(31) 


(32) 


By superposing the 3 matrix systems above, we may write 




”f • 

z y X ] 


= D 


-1 


(F . X,) (M , X,) (F , X,) 
z 1 y 1 X 1 

^F , XJ (M , X J (F , X ,) 
z 2 y 2 X 2 

(F , XJ (M . X.) (F . X J 
z 3 y . 3 X 3 


(33) 


(3 X 3) 


(3 X 3) 


(3 X 3) 
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ThuK, equatiun 16 t:an be extended tQ include three hub outputs as follows: 


i: 

V M , 
t F ( 
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M (f) 


- F ( 0 ^ 
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K.0) 

zo 

1 1 

- M (f) 
yo 1 

( = »F • »M - "f 1' 

VX^(O) 

( z y x) 
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1X3(0; 


( 34 ) 


The higher harmonic solution input is obtained by setting F , My and F to 
zero, inverting the complex transfer function matrix and multiplying: 




(X,(f)> 

1 

r-F (f)\ 

' = Ih_ , H , H 

1 F M t , 

> zo 1 

) \ 

X^(f) 

-M (f)> 
yo / 

/ 1 

X 

>. 

N 

1 X 3 ( 0 ] 


-F (o) 
f xo / 


where ’l‘ now denotes an optimal solution. 

It is seen that the transfer function generated in this estimation technique 
can be used to relate 3/ rev, 4/rev, and 5/rev harmonic blade pitching to 
similar harmonics of blade flapwise and edgewise root shears. Writing 
third, fourth, and fifth harmonic blade pitching in terms of fourth harmonic 
swashplate pitching, rolling and collective mottun in addition to solving for 
similar harmonics of blade root shears in terms of fourth harmonic hub 
forces and moments permits the direct application of equation 33. Once an 
objective function is written in terms of swashplate displacements and hub 
forces and moments, optimal swashplate inputs required to null certain hub 
vibratory forces can be derived, as in equation 35. Reference 13 presents 
details for the construction of a confidence region for the estimate of the 
frequency response function at a given frequency. In addition, special 
smoothing and filtering techniques may be necessary to improve the statis- 
tical nature of the sampled data. Reference IZ lists several frequency 
domain techniques for smoothing raw spectra. 

The anticipated control sequence for applying such an algorithm involves 
measuring 3 hub baseline vibration levels to initiate the solution procedure. 
Once the ambient or baseline levels have been determined, a single simul- 
taneous swashplate perturbation in pitch, roll and collective at 4P yields 



inlux'mation nccesdary to generate required inputs fur millin)> the bast>line 
condition. The dependency of such akobient vibration levels upon flixjht 
condition requites that the solution feathering inputs be updated in response 
to chan^e8 in rotor trim. It would be undesirable, h«jwever, to render the 
actuator system idle each time baseline levels are to be obtained for a 
solution refresh. An alternate control scheme might consist of working 
with finite differences to minimize, not the baseline response, but the 
change in hub response due to changes in rotor trim. That is, once optimal 
swashplate inputs have been derived and implemented, there exists a mini- 
mized but finite residual hub response. As these vibratory loads increase 
in response to changes in rotor trim, the control objective becomes fine of 
minimizing the "delta" sucli as to drive hub loads back to the original mini- 
mized level. This technique would, therefore, preclude idling the actuators 
each time rotor trim is changed. 

.Note here that the envisioned adaptive control network design shall not 
address alleviation of helicopter gust response. Not only arc rotors highly 
effective gust attenuators in their pres^mt form, but the low frequency spec- 
tral content of gust response quantities is not amenable to attenuation with 
fourth harmonic feathering. 

The solution algorithms presented here will bo evaluated in upcoming wind 
tunnel investigations of multiple-input, multiple-output higher harmonic 
control of NASA/ Langley Research Center. 

Conclusions 


Data obtained from a recent wind tunnel investigation of sing Ic - input, single- 
output higher harmonic control has led to the following conclusions. 

• By varying phase and amplitude of 4P collective, lateral cyclic 
or longitudinal cyclic feathering, the -4F spectral component of 
any of six rotor -transmitted hub oscillatory responses can be 
minimized for a given trim condition. 

• -IF collective inputs needed to minimize 4F hub normal forces 
induce higher peak-to-peak torsional moments and, hence, 
pitch link loads on an articulated rotor. 

• Flapwise and chordwise bending moments were fairly insensitive 
to "optimal" 4P collective inputs, on the rotor tested. 
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An invcitigation of several techniciiies for pr edictiiif; -IP swashplate inputs 
needed to minimize 4F hub vibrate; ry responses usin^ wind tunnel test data 
has generated the following con< lusionst 

• There exists an almost linear relationship between 4F hub 
responses and 4P feathering inputs. 

• Optimal single inputs can be generated from vibratory response 
data. Such inputs can be calculated from a completely general 
6-point nonlinear algorithm. However, by taking advantage of 
several key assumptions, a computationally more efficient tech- 
nique can be derived requiring only I sample response data 
points. 

• Techniques exist for treating the multiple-input, multiple-output 
mode of higher harmonic control. The effectiveness of these 
algorithms will be assessed in an upcoming wind tunnel program. 
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